Isospin coupling in time-dependent-mean-field theories and decay of isovector 

excitations 
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We show that isospin non-diagonal terms should appear in the mean field Hamiltonian when 
neutron-proton symmetry is broken. They give rise to charge mixing in the single-particle wave- 
functions. We study the Time Dependent Hartree-Fock response of a charge-exchange excitation 
which generates a charge mixing in Ca isotopes. We find an enhancement of the low energy proton 
emission in neutron-rich isotopes interpreted in terms of a charge oscillation below the barrier. 

PACS numbers: 



Atomic nuclei are a good laboratory to test quantum 
mechanics of few-body systems. The nucleon-nucleon 
(NN) interaction can be interpreted in terms of meson 
exchange and depends on the quantum numbers of the 
nucleons (spin, isospin, parity...). For instance two pro- 
tons or two neutrons cannot exchange a single charged 
pion, unlike a proton-neutron pair. The NN force can be 
divided into a central, a spin-orbit and a tensor part p| , 
the latter being mainly mediated by pions. The descrip- 
tion of the structure and the dynamic of nuclei through 
mean-field-like methods Q makes use of effective inter- 
actions such as the Gogny Q or the Skyrme Q forces, 
taking care of both the renormalization of the NN inter- 
action in the medium and of the three-body force. In such 
models, the residual tensor component is rarely explicitly 
included but the renormalization of the central and spin- 
orbit parts by the original tensor, and in particular the 
saturating effect of the latter, is incorporated on a phe- 
nomenological basis through the density dependence and 
the fitting procedure of the force. This is an important 
aspect since the tensor force contributes strongly to the 
binding and spin-orbit splittings of nuclei 5] . 

However, treated fully or partially, the exchange of 
charged pions can only be properly accounted for in a 
mean-field-like theory through the breaking of isospin 
symmetry of single-particle (s.p.) states. Indeed, includ- 
ing important parts of the residual interaction through 
breaking of symmetries is a key concept of self-consistent 
mean-field methods. However, mixing the isospin pro- 
jection of s.p. states leads to a many-body state with 
no good charge any more. Such a broken symmetry is 
to be restored in structure calculations since the nucleus 
is meant to have isospin projection (i.e. the charge) as 
a good quantum number. Sugimoto et al [f| showed re- 
cently that a Hartree-Fock (HF) 0, calculation allow- 
ing for the explicit mixing of protons and neutrons (as 
well as of parity) , and followed by the full Variation Af- 
ter Projection technique [9(, was indeed a powerful tool 
to include correlations associated with the explicit treat- 
ment of the tensor force. Such involved calculations can 
only be applied to very light nuclei however. Also, the 



complete Hartree-Fock-Bogoliubov (HFB) formalism in- 
cluding charge mixing in both the particle-hole (p-h) and 
the particle-particle (p-p) channels was derived recently 
by Perlihska et al [10] within the Density Functional The- 
ory framework. In this context, both the neutron-proton 
pairing and the correlations associated with the tensor 
force in the p-h channel could be studied. However, no 
calculation has been performed so far. 



Symmetry breaking can also affect reaction mecha- 
nisms. When correctly introduced and restored in the ini- 
tial state, it may account for the propagation of ground- 
state correlations. Moreover, reactions may induce an 
explicit symmetry breaking. For instance, the wave func- 
tion of one of the collision partners in a Charge- Exchange 
(CE) reaction becomes a superposition of isobaric nuclei. 
While the wave function describing both collision part- 
ners has to conserve charge exactly, both nuclei are en- 
tangled after the collision, i.e. the wave function obtained 
by projection on one of the two separated products have 
not a good isospin before any detection takes place. Each 
component of the isospin T3 affects the evolution of the 
others through the non linearities of the mean field. 



Mean-field methods should be generalized to take into 
account such dynamical symmetry breaking induced by 
the reaction. The aim of the present letter is to in- 
clude, for the first time, isospin-symmetry breakin g in 
Time-Dependent Hartree-Fock (TDHF) calculations 
within the Skyrme energy functional framework. As a 
first application, the effect of correlations associated with 
s.p. isospin mixing on the nucleon emission following a 
CE excitation are studied. 



In mean-field methods, allowing for s.p. isospin mixing 
translates into the appearance of non-diagonal terms in 
the HF Hamiltonian in isospin space. We derive such 
terms starting from an effective Skyrme interaction 0: 
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where 3 = 5(r x - f 2 ), R 12 = (f x + f 2 )/2, k = 

acts on the right, while its Hermitian conjugated k' acts 
on the left and P a = i(l + <Xi.<7 2 ) is the spin exchange 
operator. In Eq. Q p(r) is the total local density de- 
fined from the one-body density matrix (rsq\p\r' s'q') = 
J2 h l Phi v ' s ' ( l') L Ph{Ysq), where (fh(rsq) = (rsq\h) denotes 
the component with a spin s and isospin q of the occu- 
pied s.p. wave-function h. It is convenient to decompose 
p into its spin-isospin components acting only in coordi- 
nate space: 
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where £ 7J)Q is the antisymmetric tensor. 

The HF Hamiltonian is obtained through the func- 
tional derivative h[p] = dE[p]/dp T . One expresses the 
variation of the energy through the variation of the var- 
ious local densities defined above. One can show in 
the same way as in Ref. |l2j that their contributions 
to the HF Hamiltonian are 6p all — » GaT^i f$T a)Jf 

where the V operators act on each term sitting on their 
right, including the wave functions. Finally the HF 
Hamiltonian reads as: 
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where the a index relates to the spin degree of freedom 
and p to the isospin one (7 or rj are kept for the 3D- 
space). Whenever a (p) is 0, only the scalar (isoscalar) 
part of the density is considered, i.e. a a (f M ) is re- 
placed by the identity, <jq — l s (f = In Eq. 
densities with p = 1,2 are non-diagonal in isospin and 
will then describe the charge mixing. The HF ground- 
state energy associated with the Skyrme force defined 
in Eq. Q can be written as an integral of a local en- 
ergy density 7i(r) which is a functional of the spin- 
isospin local densities p Q ^(r) = {v\p ail \v), kinetic den- 
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and current 
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densities J 7QM (r) = i [(V 7 - V 7 ) (r|/5 QM |r')] r=r ,. In or- 
der to shorten the equations, the r dependence is not 
written explicitely in the following. The contributions to 
Ti associated with to, ti, £ 2 and £3 terms in Eq. (JTJ have 
similar structures: 

n k = c k j^ dt^r , 

a,/i,=0 

with df = bHx 2 ,df\ a ^ = df I = 2ufl, d£" | aW 40 = 1 
and, for k € {0, 1,3}, d° k ° = 3, df\ a ? = 2x k -l, 4"\i#o = 
—2x k — l and \ a ^o = ~ !■ The functions Cfe(r) are 
defined by G = |, Ci = G 2 = ||, G 3 = and 
the functional T by = = p 2 , 
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and J-^ =4p afl Ap afl ~J 7 " >J ' . The spin-orbit term is: 
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with the operators G given by Gq M = Gg M = p QM , 
G^ = 4Ap QM - G" M and 
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where the derivative 9 7 and A acts only on the func- 
tion next to them. Everything derived above is valid 
both for HF and TDHF calculations which corresponds 

to [ft|£(t)],p(t)] =ihd tP (t) [12. 

The breaking of the isospin symmetry leads to a non 
diagonal h in isospin space because of the fi )2 operators 
in Eq.Q. Through the dynamics, such a non-diagonal 
Hamiltonian will lead to an isospin oscillation analogous 
to the Larmor precession of spins in an unaligned mag- 
netic field or to particle oscillation in a non diagonal mass 
operator. In order to quantitatively illustrate the phe- 
nomenon of isospin mixing in the case of a CE excitation, 
let us use a simplified version of the complete Skyrme 
force introduced in Eq. (JTJ: 



v = [f[poo] + fAPoo]P<r) $ 
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with /[poo] = t + tfpoo and f a [poa] = t x Q + ^x 3 p^ . 
The parameters t = -1803.38 MeV.fm 3 and t 3 = 




p (fm -3 ) 



FIG. 1: EoS of pure neutron matter. Variational calculation 
with Eq. JSJ with different values of xo- 
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12912.80 MeV.fm^+W have been adjusted to repro- 
duce the empirical saturation point of infinite symmetric 
matter whereas (3 — 1/3 was chosen to have a reason- 
able compressibility Koo = 236.2 MeV. As momentum- 
dependent and spin-orbit terms are dropped, the nucleon 
effective mass m* is equal to the bare one and the scalar 
time odd current density is not generated by the force, 
while spin-orbit splittings and some magic numbers are 
not reproduced in finite nuclei. Asking for a value of the 
symmetry energy a s = 32.5 MeV at saturation density 
Psat — 0.16 fm -3 gives the relation X3 = 1.544a;o + 0.161. 
Different values of xq were tried in order to match the 
Equation of State (EoS) of infinite neutron matter pre- 
dicted by Akmal et al. through variational chain summa- 
tion methods [13j (thick-solid line on Fig. ^) . The value 
of xo = 0.2 for which the EoS of neutron matter obtained 
from microscopic calculations is well reproduced was fi- 
nally used (see Fig. In view of the commonly accepted 
value chosen the symmetry energie and the fair reproduc- 
tion of the EoS of neutron matter, one can expect the 
isovector properties of the simplified Skyrme force and 
the isospin mixing arising from it to be reasonable. In 
addition, the proton and neutron separation energies are 
close to the experimental values for the 40 Ca: S p = 8.25 
MeV and S n = 16.56 MeV. However, going far from sta- 
bility the separation energies are not as good. For 60 Ca, 
a HFB-Lipkin-Nogami calculation gives S p = 24.33 MeV 
and S n = 3.57 MeV Q whereas we get S p = 19.35 MeV 
and S n = 8.38 MeV. The full Skyrme force of Eq. [T]is 
needed to correct for this discrepancy. 

We now write the TDHF potential associated with the 
force in Eq. JSJ, focusing on even-even nuclei. In their 
ground state, only time-reversal invariant local densities 
are non-zero. Thus, s.p. states with a good spin can be 
used. This remains true during the dynamical evolution 
as long as we do not break explicitly the spin symmetry 
in the wave function. With this hypothesis, all terms 
with a a index vanish in Eq. Finally, we are left 



FIG. 2: Proton (thick solidline) and neutron (thick dashed 
line) diagonal potentials U qq and isospin coupling modulous 
\U q -q\ divided by the boost strength e at 30, 360 and 690 
fm/c (thin solid, dashed and dotted lines respectively) for 
non spin- flip CE excitations in 40 Ca (up) and 60 Ca (down). 



with the potential matrix elements: 

(f + 2f a ) + 6 9th U c 

U g - g = -\(f + 2f a ) Vq , (6) 

where U c is the Coulomb potential. The non-diagonal 
part U q - q , or "isospin coupling" is proportional to the 
isospin mixing density, rj q = ^poi+«9Po2 = J2s{ rs l\p\ rs ~ 
q) . We can now study a non spin-flip CE excitation gen- 
erated by a collective isovector boost applied at t = 
to an isospin-diagonal HF ground state \<j>o): </>(0)) = 
e~ t£T "\(f)o) where e and u denote the boost strength and 
direction in isospin space, respectively. When u is not 
aligned with the T3 axis, a CE takes place and \(f>(0)) is 
then an antisymmetrized product of s.p. wave-functions 
which mixed protons and neutrons. For N — Z nuclei the 
CE produced is quenched by the Pauli principle (and van- 
ishes if Coulomb is neglected). As explained earlier, the 
wave function after the CE excitation is a linear combi- 
nation of isobaric analog nuclei. The influence on the dy- 
namics of such a quantum entanglement in isospin space 
is incorporated in the TDHF method developped in the 
present letter. 

We have computed numerically the TDHF evolution 
after a boost an isovector applied on the ground state 
of 40 Ca and 60 Ca. The simplified Skyrme mean field Eq. 
© is considered. lie « 1, the diagonal potentials (U qq ) 
are almost constant in time while the isospin coupling is 
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FIG. 3: Ratio of the numbers of emitted proton following 
an isovector excitation by the n and T3 operators as func- 
tion of the proton energy E divided by the barrier B for the 
40-48-54-eo Ca ( solidj d as hed, dotted, dotted-dahed lines re- 
spectively). 



linear in e. The resulting potentials at different times are 
plotted in Fig. |21 The small variation in time of |£7 q _ q | 
comes from the difference in energy of the excited modes. 
|t/ g _ g | is more intense in 60 Ca than in 40 Ca as expected 
from the quenching of CE processes in the N=Z nuclei. 
The isospin coupling is also more peaked at the surface 
in 60 Ca because of the isospin mixing in the 2p-\f shells 
and decreases exponentially at the barrier. 

Let us now present a possible effect of the isospin cou- 
pling on the low energy proton emission. The particle 
emission looked at in the present case originates from 
the entangled state prior to any detection of the reac- 
tion partner projection on good isospin). An emit- 
ted proton (neutron) can come from an incident pro- 
ton (neutron) or an incident neutron (proton) which 
has exchanged its isospin via the isospin coupling. At 
low energy, where the proton emission is hindered by 
the Coulomb barrier, one can expect an enhancement of 
the proton emission from neutron which exchange their 
isospin at the surface. This is possible because the isospin 
coupling does not fully vanish at the surface (see Fig. [2Jl . 
The excitation studied above does not usually populate 
states above the nucleon threshold. To reach such high 
excitation energy we applied on the HF ground state an 
isovector-monopole boost e~ ler Tu . Due to the r oper- 



ator, an important part of the strength goes to the 2Hlu 
non spin-flip Isovector Giant Monopole Resonance. Let 
us compare the numbers of protons emitted at an energy 
E after a boost generated by r 3 (noted Np) and by t\ 
(ATi). Fig. shows the evolution of N^/N^ - 1 obtained 
for 40-48-54-60 C£L We choose £ = o.Ol fin" 2 which is 
small enough to be in the linear regime. 

For a T3-boost, there is no isospin coupling whereas a 
ri excitation produces an isospin coupling slightly more 
peaked at the surface than in Fig. [5] due to the r 2 term 
in the boost. In this case the expectation value of the 
charge is changed by A{Z) = -0.0044 in 40 Ca and 0.064 
in 60 Ca. Fig. |21 shows a strong enhancement of the pro- 
ton emission at energies well below the barrier except for 
the 40 Ca for which the isospin coupling remains small. 
Around the barrier this effect is hindered by the differ- 
ence in the excitation energy spectra associated with the 
inelastic and the charge exchange excitations. This en- 
hancement of the low energy proton is expected to be 
a general consequence of isospin coupling. In fact we 
expect this effect to be strong in all nuclei with large 
proton/neutron asymmetry. 

To conclude, we have shown that isospin coupling ap- 
pears in the HF mean field when proton-neutron sym- 
metry is broken. Its expression was derived for a full 
Skyrme-like force. We used a simplified Skyrme force ad- 
justed on the isospin properties of nuclear matter to get 
an approximation of the isospin mixing. As a first appli- 
cation we performed TDHF calculations using this force 
to study the effect of charge-exchange excitations on nu- 
cleon emission. We found an enhancement of the proton 
emission below the barrier for the neutron rich Ca inter- 
preted in terms of nucleon charge oscillation. More quan- 
titative predictions are needed with a 3D-TDHF code 
with a full Skyrme force to study complete reactions, i.e. 
not only the evolution of one nucleus wave function. The 
inclusion of the spin degree of freedom would also allow 
studies of Gamow- Teller transitions and isospin and spin- 
flip giant resonances. For the latter case however and for 
any nuclear structure study in general, it is necessary to 
keep the good value of the charge in average during the 
time evolution and to project on it after. 
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